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Abstract 

The wave functions corresponding to the zero energy eigenvalue of a one- 
dimensional quantum chain Hamiltonian can be written in a simple way us- 
ing quadratic algebras. Hamiltonians describing stochastic processes have 
stationary states given by such wave functions and various quadratic alge- 
bras were found and applied to several diffusions processes. We show that 
similar methods can also be applied for equilibrium processes. As an exam- 
ple, for a class of g-deformed O(N) symmetric antiferromagnetic quantum 
chains, we give the zero energy wave functions for periodic boundary condi- 
tions corresponding to momenta zero and n. We also consider free and various 
non-diagonal boundary conditions and give the corresponding wave functions. 
All correlation lengths are derived. 
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I. INTRODUCTION 



Quadratic algebras and their representations have being extensively used recently in 
order to study the probability distributions of steady states of one-dimensional stochastic 
processes with open boundaries or on a ring . The basic idea is that if the Hamiltonian 
of a quantum chain which gives the time evolution of the system, has eigenvalue zero, the 
ket wave functions which are related to the steady states probability distributions have a 
simple expression in terms of a certain quadratic algebra determined by the bulk rates. This 
algebra has representations fixed by the boundary conditions, the corresponding matrices 
act in an auxiliary vector space. All correlation functions can be computed from these ket 
wave functions. The aim of this paper is to "import" these techniques to equilibrium statis- 
tical physics and stress the limitations and differences. For stochastic processes the lowest 
eigenvalue of the Hamiltonian which gives the time evolution of the system is zero. This 
is not the case for most of the Hamiltonians which are interesting in equilibrium problems. 
Therefore the possible applications of the algebraic approach to ground-states is bound to be 
more limited. Another difference is that in equilibrium and periodic boundary conditions, 
the ground-state can have momentum non-zero (is not translational invariant). This can't 
be the case for stochastic processes since the components of the ground-state ket vector 
have to be positive numbers (they are probabilities). Another difference appears when we 
want to calculate correlation functions which are expressed in terms of vacuum expectation 
values (implying the bra AND ket vacua). As we are going to see the expressions of the 
correlation functions are very similar in both cases. Actually the quadratic algebras ap- 
proach was implicit already used in equilibrium problems where it is known as the matrix 
product approach [Q-fJ. The matrices used are in fact representations of certain algebras. 
We hope to convince the reader that the algebraic approach is not only more aesthetic but 
more powerful since it makes contact with known results obtained in mathematics. Finally, 
we would like to mention that matrix product approach has been used as an alternative 
to the density matrix renormalization group method [0|§. How the methods presented in 
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this paper can be applied to this problem is an open question. The application of quadratic 
algebras to zero-energy states is presented in Sec. 2. Much of the contend of this Section is 
already known. What is new is how to handle zone boundary states which have momentum 
Ti. In Sec. 3 we give an application. The idea is simple: in the study of quantum groups ||, 
in order to find the non-commutative manifold in which they act, Reshetikhin et al have 
introduced projector operators out of which one can build quantum chains having the quan- 
tum algebra as symmetry. Moreover, one gets for free a quadratic algebra (the manifold of 
the quantum group) which can be used to write the zero energy eigenfunctions of quantum 
chains build using the projector operators. These chains are not exactly integrable in the 
Yang-Baxter sense. We have considered, as an example, the O(N) case for which we get 
an iV-state Hamiltonian. The quadratic algebra turns out to be the g-deformed Clifford 
algebra. In the special case N = 3 and q = 1 one recovers the model with valence bond 
ground state (VBS) of Affleck et al 0. (The q-deformed case can be found in Refs. fTOB). 
The N = 4 case is discussed in the Appendix A, it is a special case of the extended Hubbard 



model O] . The Hamiltonians we consider can be mapped into quantum spin ladder models 
| T2| and find applications in this context. We are going to show that for periodic boundary 
conditions and an even number of sites, we find a unique momentum zero ground-state. For 
N even, we also find one zone boundary state. For free boundary conditions, we find 2 N ~~ 1 
ground-states. This degeneracy can be lifted adding boundary fields. In Sec. 4 we show how 
to choose the boundary conditions in order to get an unique vacuum. The boundary terms 
break the symmetry of the quantum chain. The calculation of all the correlation lengths (for 
any N) is presented in Sec. 5. It is shown that for large N the correlation lengths diverge. 
In Appendix B we show how to compute the correlation function for some parity violating 
operators, appearing in the case where N is even. This problem is interesting in the case 
of periodic boundary conditions when the ground-state is twice degenerate even for a finite 
number of sites. The conclusions can be found in Sec. 6. 



3 



II. ZERO ENERGY STATES AND QUADRATIC ALGEBRAS. 



The application of quadratic algebras to the zero energy ket wave functions for diffusion- 
reaction processes is well known [ffj,0,H, in this section we will do a trivial extension to 
equilibrium processes and show how to compute correlation functions. We consider a most 
general one-dimensional quantum chain with N states, L sites and nearest-neighbour two 
body interactions. The Hamiltonian is: 

L-l 

H = Y,H k + c + n. (i) 

k=l 

The bulk terms {k — 1, . . . , L — 1) and the left and right boundary terms are: 



N 

E 



H k = E r#3Z°<i (2) 



N N 

£= J2 tfEt, K = £ R%E{\ (3) 

a,/3=l a,(3=l 

Here Et basis for N x N matrices on the k-th site: 

= 5 ai 5 p5 (a, p, 7 , 5 = 1, . . . , N). (4) 
We will assume that H has at least one eigenstate of energy zero 

H\0>=0 <0\H = 0. (5) 

Our aim is to describe the bra < 0| and ket |0 > states in a simple way. In order to do 
that, we consider two associative algebras defined by the bulk interaction: 

N 

T^x a Xfs = x^X s - X^x s (6) 

a,/3=l 
N 

E T "sy^ys = VcXp - y a yp- (7) 

7,5=1 

If the bulk part of the Hamiltonian is not symmetric, the two algebras are different. Each 
algebra has 2N generators x a , X a and y a , Y a ,( a = 1, . . . , N), respectively. We define two 
Fock-like representations of the two algebras: 



N N 

< V K \(X a - J2 L^Xfs) = (X a + R&pWk >= (8) 

0=1 (3=1 
N N 

< v B \{Y P - £ Liva) = (Yp + Y, Riy a )\w B >= 0. (9) 

Here < Vk\, \ Wk >,< Vb\ and \Wb > are the bra and ket reference states defined by the 
equations (§) and in AUXILIARY spaces. We make now the connexion between the two 
algebras and the zero energy eigenstates of the Hamiltonian. The basis in the ket vector 
space in which the Hamiltonian acts is: 

u ai u a2 ...u aL (at = l,2,...,N) (10) 

the iV-dimensional vector u ak is in the k-th site and has the component equal to one and 
the others zero: 

(««J/j = *a^ (/? = 1, 2,..., TV). (11) 
We denote the basis in the bra vector space in which the Hamiltonian acts by 

«•••<■ (12) 

The scalar product is obviously 

< u l k u P k >= ( 13 ) 

One can prove |T3 that the unnormalized bra and ket vacua can be written using the two 



quadratic algebras: 



N 



|0 >= Yl <V K \x ai ...x aL \W K >u ai ...u aL (14) 



ai,...,a L =l 

N 

<0|= £ <V B \y ai ...y aL \W B > u T ai ...ul L . (15) 

otl,— ,ai,=l 

Notice that the generators X a and Y a don't appear in the expressions of the wave functions. 
One can also show that the quadratic algebras exist, and that one can find representations 



satisfying the conditions fl8|) and ([|). Moreover, one can show that all the zero energy wave 



functions can be obtained in this way [ II] . In the case of periodic boundary conditions, and 
translationally invariant zero energy eigenfunctions, one can use the expressions (fH| ) and 
(|T5|) making the substitution: 

<V K \...\W K >^Tr(...) <V B \...\W B >^Tr(...) (16) 

provided that the algebra has a trace operation. 

As opposed to the case of the Hamiltonian with open boundaries, for periodic bound- 
ary conditions, it is not clear in which cases one obtains in this way all the zero energy 



eigenfunctions. A simple counter-example was given in Ref. [15] in which it is shown that 



there are zero energy eigenfunctions which can't be obtained using the algebraic method 



given by equation (|T6|). On the other hand, examples are known [16j where indeed all the 
eigenfunctions are obtained. 

Ground-state wave functions can correspond to zone boundary states (momentum 7r). 
One can show that if the algebra ([|) has the Str operation with the properties: 

Str(x Ql x Q2 • • • ■J-qi) Str(x Q , L x Q , 1 x Q , 2 • • • ^o^—i) 

Str(A Q , 1 x Q2 • • • ^'a i, ) Str(x Qi X ai x a 2 . . .x aL1 ^j (1^) 
than the ket vector 

L 

|0 >= Str(x Ql . . .x aL )u ai . . .u aL (18) 

ai,...,a L =l 

satisfies equation @ and it is obviously a zone boundary state. Similar expressions can be 
used for the algebra (0) and the bra eigenvector. The Str (called supertrace) operation is 
taken from the theory of superalgebras and it implies that the x a and X a are odd generators 
in this algebra. In particular if in the algebra @ one takes X a c-numbers (this is often done 
for diffusion processes [[|), the algebra can't have the Str operation. In Sec. 3 we will show 
in examples how the Str operation works. As for translationally invariant ground-states it 
is not known if all of the zone boundary states can be obtained using equation §Li 
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Before showing how to compute correlation functions, let us see what are the conse- 
quences for the quadratic algebras of the existence of a symmetry of the Hamiltonian. Let 
us assume that the operator: 

L N 

A = E E A^EfT (19) 

k=l fi,v=l 

commutes with the bulk part of the Hamiltonian, i. e., 

[A, E H k \ = 0. (20) 

k=l 

Simple arithmetics gives the relations: 

TV 
«,/3=l 

(A^x^)X 5 + x y (AsnXn) - (AypX^xg - X y (AspXp). (21) 

This relation gives a set of simplified algebraic relations among the generators of the algebra 
and at the same time, shows that the generators are tensor operators. (A relation similar to 



(pl|) can be obtained for the generators y a and Y a ). As an example, let us choose A n = 1 



and all the other matrix elements zero in ([L9|) . Using ( pT|) one obtains 



N 



E(r^xix a + r$x aXl ) = 6 lA ( Xl X 5 - X x x s ) + 5si(x J X 1 - X jXl ). (22) 



a=l 



Similar relations can be obtained in the case of quantum algebra symmetries when the 
operator A has not the simple expression (^9[). We now show how to compute a two-point 
function. This calculation is interesting when the ground-state energy is zero. Consider two 
local operators P r and Q s on the r and s sites. They act on the basis fllCf ) as follows: 

N N 
PrU ar = E P Pr,a r Up r ; Q s U as = E Qp s ,cc s Up s . (23) 

/3r = l P,=l 

We want to compute the expression: 

Gr ^<m$^ (24) 

where < 0| and |0 > are given by equations (|HD and (p~5| ) and Z is a normalization factor 
coming from the fact that (0) and flTSP give unnormalized wave functions. It is useful to 



define the following quantities (all related to the auxiliary space) 



N 

c = Y, x *®y* ( 25 ) 

a=l 

N N 

P = P aP x P ®ya, Q = Qa,f3X(3 ® y a (26) 

0,(9=1 a,(9=l 

and 

< Vb|®< Vkr| =< V|; |W >= > ®\W B > . (27) 
Using equations (G5J)-(B7[), the two-point function (E^) has the following simple expression: 

G r , s = - < V\C r - l PC s - r - l QC L - s \W > (28) 

Zj 

where 

Z=<V|C L |PU>. (29) 

Notice that C plays the role of a space evolution operator in the auxiliary space but the 
analogy with a quantum mechanical problem can't be pushed further since < V\ and \W > 
are not eigenfunctions of C. Nevertheless one can see that the spectrum of C gives all 
the correlation lengths. For periodic boundary conditions, we have to make the following 
substitution: 

< V\...\W >^Tr(. ..) 

< V\ ... \W >-> Str(. ..) (30) 

for translationally invariant states, or for zone boundary states, in equations (|28|) and (p9|). 
Let us observe that the expressions fl2~5])-(29"D are similar to the ones one obtain for stochastic 



processes . The difference is that instead of dealing with only one algebra (given by 
equation (|)), one has the tensor product of two algebras. If the algebra (0) has a one- 
dimensional representation, (this is always the case for diffusion processes with exclusion for 
example |IJ), the correlation functions computed using the ket vector only or the bra and 
ket vector (vacuum expectation values), coincide. Expressions like (^5|)-(p9|) have been used 



in a different context in the matrix product approach to the density matrix renormalization 
group method |7|]. In this case the x a are matrices obtained using the variational method 
and not using quadratic algebras defined by the Hamiltonian using equation (|6|). Besides 
they have to satisfy the condition: 

JV 

53x a a£ = l. ( 31 ) 

a=l 

As we are going to see in the next section, this condition is not necessarily fulfilled in our 
applications. 

III. Q-DEFORMED O(N) SYMMETRIC, iV-STATE QUANTUM CHAINS. 

The quantum chains describing stochastic processes are given by non-hermitian Hamilto- 



nians which always have zero as lowest eigenvalue. The quadratic algebra always exists fI3 
and the problem is to find representations of the algebra. In equilibrium problems one is in- 
terested in hermitian Hamiltonians which in general don't have zero as the lowest eigenvalue 
and therefore one has to find Hamiltonians which have this property. In order to illustrate 
the method, in this paper we have chosen an easy way: using known results in the theory 
of quantum groups. In this way we get not only hermitian quantum chains which have zero 
for the ground state energy but also quadratic algebras with known representations. 

A. The bulk Hamiltonian. 

Reading the paper of Reshetikhin et al one can notice that there are several expressions 
of the form (|6|) with the X as equal to zero. We will choose the one where are projector 
operators of rank N(N + l)/2 — 1 for the g-deformed B(n) series (N — 2n+ 1) and D(n) 
series (N = 2n). The x a are the generators of the non-commutative algebra of the manifold 
where the quantum groups act. Similar expressions for the Sp(n) and Osp(m/n) algebras 
and superalgebras can also be obtained (T7J. In the present paper, we confine ourselves to 
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the q-deformed O(N) case. As we will show we will use these projectors in order to write 
Hamiltonians for quantum chains. The projector operators have the following expressions: 



P 



(+) 



N 



1 



^aj3 77,70 t-,5/3 _ 

k ~ L 75 ^k+i - - 1 

a,/3,7,5=i y ' y 



? E « + - g" 1 ) E £f 



N+l N+l N+l N+l 
' 2 

k+1 



N 



+0N,2n+lE k 2 2 E, 2 



t 1 E *T + E £f < a 

a,/3=l a^/3,/3' 



-2L N 

£ i?r'< a i -jwr E - - g- 1 ) E K ,/3 <'iy 

L 2 J<? «,/3=l 



(32) 



a^a' L 2 J<? «,/3=l a>/3 

where g is a deformation parameter (taken real in this paper) and we use the notation 



n 



q n — q v 



q-q 

In equation (|32|) we also denote 



— and a' = N+l -a (a = l,...,N). 



(pi, Piv) = (n - -, n 0, -n + -), 



for iV = 2n + 1, and 



(33) 



(34) 



(Pi, • • • , Pn) = (n - 1, n - 2, . . . , 1, 0, 0, -1, . . . , -n + 1) 



(35) 



for N=2n. By definition we have 



(P; 



(+)\2 _ p(+) 



(36) 



Since the matrix in (|32|) is symmetric, i. e., = T^, the two associated algebras to 



the projector (|32]): 



JV JV 

E F^x a xp = and ^ T^y a y^ = 

a,/3=l «,/?=! 



(37) 



are identical and therefore we give only one of them. It is convenient to denote (for obvious 
reasons); for N=2n: 

%1 = 0"ni %2 = 0"n—li • • • i %n = = O'li %n+2 = a 2 J • • • ) X 2n = a n (38) 

and for, N=2n+1 
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X\ O n , %2 O n _i, . . . , X n Oj, Xn+1 



Vs + s 



V - + H 

= Zj, X n +2 — "i , • • • , X2n+1 — a Tl 



(39) 



where s = +/q. The 2n g-deformed fermionic creation and annihilation operators a a , a+ and 
the "7 5 "-type generator E satisfy the following relations: 

qapa a + a a a/3 = (/3 > a) 
qa a^ + a+ap = ((3 > a) 
Sa a + ga a S = 0, S + = S 
a a a+ + a+a a = qa a+1 a+ +1 + g -1 a+ +1 a Q+1 
gaiety + g a^ai = L . 



1 < a < n- 1) 



(40) 



The above algebra has a central element : 

and an obvious representation is: 

a fc = l<g>l<g>---<g>a<g> s CT V 2 ® s° z a z <g> ■ • 



(41) 



s° z a z , (k — 1, . . . ,n) 



Y, = s a a z (g) s a a z (g) ■ ■ ■ (g) s a a z , 



(42) 



with 









f 0^ 






"I 


a = 




, a + = 












^0 oy 




I 1 °y 




^° 


- 1 / 



(43) 



In the first line of (|42|) the operator a is in the fc-th position, and the operator s u a z appears 
in the positions k + 1, . . . , n. The fact that the algebra has finite-dimensional representations, 
makes all calculations much simpler (see Sees. 5 and 6) as compared with the cases when the 
algebra has infinite dimensional representations. Notice also that for q ^ 1, the generators 
x a do not satisfy the relation ([IT]). We make now the connexion between the projectors (|3"2"D 
and the quantum chain ([TJ). Since the lowest eigenvalue E of a projector operator is zero, 
we can choose in ([!]): 
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H k = Pt +) (44) 

Notice that is hermitian and therefore since as mentioned its lowest eigenvalue is zero, H 
for periodic or free boundary conditions has also zero as its lowest eigenvalue. The problem 
of other boundary conditions is going to be discussed in Sec. 4. 

B. Ground-states for periodic and free boundary conditions. 

We start with periodic boundary conditions. We first consider zero momentum states. 
Using equation ( |T4| ) (together with the substitution given by (|T6| ) as well as the represen- 
tation (f4*2"D we get for all N, one single ket vector of energy zero for L even and none for L 
odd. This result is confirmed by the spectra obtained from the numerical diagonalization 
of several Hamiltonians (various L and N). This check was necessary since as mentioned 
in Sec. 2 there is no theorem which assures us that there are no zero energy eigenf unctions 
which are not obtained using the algebraic procedure. We now look for zone boundary states 
and therefore look for a definition of the Str operation such that the relations (|i~7D and (|i~8l) 
are satisfied. We consider the matrix J defined by 

J = cr z ®a z <g> ■ ■ ■ ®a z . (45) 

A vector 

\v >= \vi > ®\v 2 > ® • • • ® \v n > (46) 

is called even (odd) if it is an eigenvector of J corresponding the the eigenvalue +1 respec- 
tively -1. A matrix is called even if it takes an even (odd) vector into an even (odd) vector. A 
matrix is called odd if it takes an even (odd) vector into an odd (even) vector. For example, 
the matrices a>k in equation ( [42]) are odd but the matrix S is even. Consider now the matrix 

A = A l <g> A 2 ® ■ ■ ■ <g> An. (47) 

We define 
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Str(M) = Tr(JM). (48) 

It is easy to check that if A and B are odd matrices, than 

Str(AB) = -Str(SA). (49) 

If one of the two matrices is even and the other one is odd 

Str(AB) = 0. (50) 

If the two matrices A and B are even 

Str(AB) = Str(BA). (51) 

From this properties we learn that in order to satisfy the relations ( |P7|) (the relations in 
the second line of ( 1~8|) are automatically satisfied since X a = 0), the x Q 's have to be all 



odd generators. This excludes the case of N = In + 1 because of the appearance of the 
sigma generator which is even. For L odd and N = 2n all the supertraces are zero and 
again we can't obtain a boundary state which is physically correct. For N and L even 
we expect therefore a unique zone boundary state. This is what is also seen in numerical 
diagonalizations for all L except for q = 1 and small values of L where something subtle 
happens. We illustrate the phenomenon taking N = 4. Using equations (|4"2"D , fl4"5]) and (0Sp 
we obtain: 

Str(a^a 2 ) = 0, Str(a^ai) = q — q^ 1 , (52) 

which would imply that for q — 1 and L = 2 there are no zone boundary states. Actually 
there are two of them which can be obtained taking instead of J given by equation (]45|). 
two alternative expressions: 

l®a z or a z ®\. (53) 

These expressions can't be used however for monomials with more than two generators (the 
property ( |i~7[ ) is not valid anymore). In the spectra for periodic boundary conditions as seen 
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in numerical diagonalizations, there are no zero energy states besides those mentioned above. 
The existence for N even of a degenerate ground state, one of positive parity (momentum 
zero, obtained with the help of the Tr operation) that we denote by |0, + > and one of 
negative parity (momentum 7T, obtained with the help of the Str operation that we denote by 
|0, — >), allows for the existence of correlation functions of operators which break parity. For 
example one of the operators P or Q in equation (|25| ) can break parity. In the Appendix B we 
show how to compute the correlation functions for this case (one considers matrix elements < 
0, — | • • • |0, + > for example). A somehow similar problem occurs in spontaneously dimerized 



spin ladders |]18| . We would like to stress that in our case the degeneracy of the vacuum 
takes place even for finite number of sites. 

We now consider free boundary conditions. An inspection of equation (H) shows that it 



brings no constrains therefore instead of equation (|T4"D we have: 

N 

|0 > ^ ] %ati ■ ■ ■ •^aL^ai ■ ■ ■ ^a^i (^4) 

where the various independent monomials (words) in the algebra are regarded as a basis 
in a vector space. Each component of |0 > in this basis gives a zero energy eigenfunction. 
Therefore for both L even and odd we get 2 1 states. This result was obtained counting 
the independent words. For small values of L, the degeneracy can be smaller since higher 
degree monomials might not yet have appeared. For example for N = 4, and L = 2 the 
degeneracy is 7 instead of 8 but for L = 3 one obtains already 8. 



IV. BOUNDARY CONDITIONS COMPATIBLE WITH THE QUADRATIC 

ALGEBRAS 

The boundary matrices £ and TZ (we will choose them hermitian) have not only to be 
compatible with the quadratic algebra (see below), but have also to leave the value zero as 
the lowest eigenvalue. This property is warrantied if the lowest eigenvalues El and Er are 
also zero. This follows from the relation: 
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E H >E L + (L- l)E + E R . (55) 

where Eh and E are the lowest eigenvalues of H and H^. 

Since for the g-deformed O(N) symmetric quantum chains defined in the last section, 
the algebras (|5|) and (0) with X Q and Yp equal to zero are identical, we have to find the 
matrices £ and 1Z as well as the vacua of the auxiliary spaces such that the following relations 
(obtained from equations ©-©) are satisfied: 

N N 

^2R^x p \W K >=0; ^2RPx a \W B >=Q, (56) 

(3=1 a=l 
N N 

y £L^<V K \x p = 0, "£L a p<V B \x a = 0. (57) 

(3=1 a=l 

We have taken the same representation for the two sets of Clifford generators x a and y a . 
We now show that the solutions of equations (|5"Tj) can be obtained from those of equations 
(|56|). We take the transpose of the two equations d57|): 

N N 

]TL^J|V£>=0, Y. L Pl\ V B>=®, (58) 

13=1 a=l 

and since from equations (|38|)-(|39D, we have = x a < = X]y + i^ a , we can rewrite the equations 
(|58l) as follows: 

N N 

Y t L%x l3 \V^>=0, Y, L i x «\VB >=0. (59) 

13=1 a=l 

We can compare now the equations (|56| ) and ( |59D and deduce that for any solution R@ of 
(p)Bj) (there are many of them) one gets a solution for L" a : 

= <• (60) 

One can use of course one solution of equations ([56]) for R 1 ^ and another solution to get 
using equation (|60D . One can easily show that the solutions of (|56|) have a factorized form: 

R a 3 = rejp. (61) 
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It is convenient to choose the following basis in the auxiliary vector-spaces in which x a and 
y a (replaced formally by x a ) act (see equations (BH) and (f42|)-(f4"3"|)): 



\W K >-- 



L 

n 



\W B >= (f[ 



i=i Jl + 77,. 



p;2 



v 1 ; 



v 1 ; 



<8> 



V 1 / 



<g> • • • ® 



V 1 / 



'In 

V 1 / 



(, \ 

'In 



(62) 



(63) 



(0)-® 

We are going to consider separately the cases N — 2, 3, . . . , 6 in order to illustrate the 
structure of the solutions. As we are going to show for iV = 2, the values of 771 and fji 
are fixed and besides a common factor, the matrix elements of 1Z contain no parameters. 
For iV = 3 and 4, the parameters 771,772 respectively 771,772 are free and 1Z is given by the 



parameters of the wave functions (32)-(|63|) and a common factor. For N = 5 and 6, a new 
phenomenon appears. The wave function (|62|) is given by the free parameters 771,772,773, and 
the wave function (|6"3"| ) is specified by the corresponding parameters 771,772,773. 71 depends 
now not only on the parameters of the wave functions but on supplementary free parameters. 
This implies that different boundary conditions are compatible with the same wave functions 
(0)-(HD- We now consider the boundary conditions for some values of N. 
N=2 



Since this is a very simple (and trivial) case, we discuss it in detail. From equations fl32|) 
and (|44]) we get: 



H k = -{alal +l + l). 



(64) 



This implies that for free boundary conditions the ground-state is twice degenerate, with 
antiferromagnetic ordering. This degeneracy is a consequence of the existence of two inde- 
pendent words in the 0(2) algebra: Xix 2 ■ ■ ■ x±x 2 and x 2 Xi ■ ■ ■ x 2 x±. Demanding that 1Z is 
diagonalizable, equations Q56p have two solutions: 
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and 



where 









(A 




(A 


11 = r A; | 


w K 


>= 




• 


\w„>= 










\ / 




i ; 
\ / 












f°) 


TZ = rB; | 


w K 


>= 






\w B >= 














k 1 j 
\ / 




(l 


°) 




'o o\ 




A = 






, B = 








v° 






v° V 





(65) 



(66) 



(67) 



and r > 0. 

We consider separately the two cases: 

a) 1Z = rA, £ = I A with r, I > 0. For lattice size L even, we obtain no zero energy 
eigenstate. The matrix elements (]14])-(|T5|) vanish. For L odd, one obtains an unique zero 
energy ground-state. 

b) 1Z = rA, C = IB For L even, one obtains an unique zero energy ground-state and 
none for L odd. 

N=3 

The solutions are : 



f\ = ST}1, f 2 = -Vs + S- 1 , f 3 = -(S77i) 1 



e x = sf}i, e 2 = -Vs + s r , e 3 = -(sfji) 1 



(68) 



and the eigenvalues of TZ are: zero two times and 

1 



1 + 



s + s~ 



(69) 



this relation imposes r > 0. Since we want 1Z symmetric, we take (see equation (51)) r\\ = fji 
and therefore: 



f a (« = l,2,...,iV). 



(70) 
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N=4 

One gets: 

fi = srfr, f 2 = -?7i, h = -Vi 1 : U = ( 71 ) 

The matrix 1Z has three eigenvalues zero and one equal to: 

r^y + vf + Vi' + is^)- 2 ]. (72) 

Notice that for N = 3 and 4, the parameters of the vacua and r determine the 7Z matrix. 

This is abound to change for larger values of N. 

N=5 

The solution is 

/ s _ s -i\ 2a 
/1 = T) 2 s 1 + a —— — ), f 2 = -77i(l - a), /; 



s + s V v s + s 



1 

1\ 



/ 4 = -^(a + 1), / B = V 1 [ 1 + *—^zi) , (73) 
where a ia an additional free parameter. 1Z has now four eigenvalues zero and one equal to 

N=6 

One gets: 

fi = V3S (a + S * j , / 2 = -r/ 2 - j , / 3 = 

/4 = -^, /s = -ife 1 (a + ^y^) , h = -s- V (« + ^y^) (74) 

with a arbitrary. 7£ has now five eigenvalues zero and one equal to rJ^Li ■ 

Notice that for iV = 5 and 6 the /'s depend not only on 771,772 and 773 but also on the 
supplementary parameter a. This implies that the same wave function can be used for 
different boundary matrices. We also notice that taking r positive makes sure that the 
lowest eigenvalue stays zero. One can obtain 1Z matrices with only non-vanishing element 
on the diagonal (like in equation fl67|)) taking one of the rji equal to zero or infinity. This 



remains valid for any N. For larger values of N the number of free parameters increases 
and it is certainly not our purpose to give here the general solution. We would like to stress 
that for N > 2, the boundary conditions can break all the symmetries of the Hamiltonian. 



V. DIAGONALIZATION OF THE C MATRIX AND CALCULATION OF THE 
CORRELATION LENGTHS OF THE Q-DEFORMED O(N) QUANTUM CHAIN. 



It is necessary to have a new look at the expression (|28 ) of the two-point correlation 
function. In the last section we have shown how to get the bra and ket vacua (< V\ and 
\W >, respectively) in the auxiliary spaces. In order to proceed further, one has to find the 
similarity transformation S which diagonalizes the matrix C: 



c = sc D s-\ 



(75) 



The matrix C is given by equation (|25|). In this equation the x a and y a are the generators 
of the two identical algebras (see (|38|)-(P0D) having the representation (p^)-(|43|). We will 
consider separately the cases where N is even or odd. 
a) N=2n 

It is convenient to write C as a four-state Hamiltonian with n sites in the auxiliary space: 



C {N) = E 1 F 2 F 3 ...F n + E 2 F 3 F 4 . . . F n + . . . + E n 



where the matrices Ei and Fi act on the i-th site and have the expression: 



(76) 



E = a®a 



(1,1) (1,2) (2,1) (2,2) 

/ 1 ^ 





10 



(77) 
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F = s a a z ®s a a 2 



\ 



(1,1) (1,2) (2,1) (2,2) 

q 
0-100 

0-10 

q- 1 



(78) 



where 



l] 



a" 



1 
-1 



(79) 



The basis vectors in the tensor products (|77|)-(|78D correspond to the two-dimensional repre- 
sentations used in In this basis the vacuum \W > in the auxiliary space has the 
expression: 



\w >= v {1) ® v {2) <g> v {n) 



(80) 



where 



V 



Notice the recurrence relation: 



(l + ?? 2 )(l + ?f) 











v 1 ; 




v 1 / 



(81) 



c{ N + 2) = C W F + ^ 



+ 1 



(82) 



that we are going to use later on. If q = 1, the diagonalization of C is trivial since E and F 
commute. Using the similarity transformation 



U 



V2 



' 1 -1 ^ 

v 7 ^ 



i 1 1 



(83) 
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we have 

E D = UEU~ l = \{° z ® 1 + 1® a 2 ); F D = U- 1 FU = a z ®a z . (84) 
It is convenient to write Cd as a one-dimensional two-state spin chain with 2n sites: 

C { D ] = + <*S)(0M) • • ■ «-!<) + 

(oj + a|)(a 5 V 2 ) ■ • • «_!<) + • • • + «_ x + <)]• (85) 

In order to simplify the expression (j55[), it is useful to look at C as a function defined on 
the abelian group Z 2 ® Z 2 ® ■ ■ ■ ® Z 2 = (Z 2 )® 2n . In order to do so, we write: 

o> k = (-l)« (ft = 0,1). (86) 
Using this notation, instead of equation fl85|) we obtain: 

cf) = + (-i) e2 )((-i) e3 (-i) e4 ) • • • ((-l)— (-i) e2 ") + 

((_!)* + (-l)-)((-l)^(-l)^) ■ ■ ■ ((-!)— (-!)«-) 

+ + + (_!)«-)]. ( 87 ) 

We make now the change of variables: 

u;i = ei + (e 3 + e 4 ) H h (e 2 „_i + e 2n ) 

^2 = £2 + (e 3 + e 4 ) H h (e 2 „_i + e 2n ) 

^3 = e 3 + (e 5 + e 6 ) H h (e 2n _i + e 2n ) 

u)a = e 4 + (e 5 + e 6 ) H h (e 2n _i + e 2n ) 



w 2n -\ — e 2n _i 

^2n = C2n- (88) 

Notice the identity: 

(_!)£<=! = (_i)E<=^ (89) 
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that we are going to use shortly. With the change of variables fl88|), instead of the expression 



37|), we get: 



N 



N 



c 



8=1 1=1 



(90) 



From equation ( P0| ) it results that the spectrum of CW for N = 2n is the same as that of 
the ^-component of the total spin S z for 2n spins | . Therefore the eigenvalues are n — m 
(m = 0, 1, . . . , N) with a degeneracy given the binomial coefficient C^~ m . 

We now consider the case q ^ 1. We are going to use the recurrence relation fl76|). We 
first make a change of basis (see equations ( ff7|) -( [f8|) ) in the four-state chain with n sites: 



(2,2) 



;i,2) 



(2,1) 



(91) 



and denote by u\ , the basis vector on the A;-th site having the z-th (i = 1, 2, 3, 4) component 
equal to one and the others zero. In this basis E and F have the expressions: 



E 



^0 10 0^ 



10 



In the same basis, for N = 2, one has 



V 



/ 



1 q ^ 

o?- 1 

0-10 
yO 0-1 



(92) 



(93) 



and the eigenvalues (eigenfunctions) are: 



V2 



-1, 



(1) (2) 



V2 



; o, 



Assume that ^1 written in the basis 



; o, 



4 1] 



(94) 



(^ = 1,2,3,4), 



is an eigenfunction of corresponding to the eigenvalue A. We now consider the four 
wave functions 
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*fM n+1) (< = 1,2,3,4), (95) 
and act with C^ N+2 ^ on them using the relation (|82"D . We obtain: 



(JV+2WJV) (n+1) _ A^(iV) (n+1) 
(AT+2WA0 (n+1) _ A , T ,(A0 ,(n+l) 



Cr +Zj ^ v; < ; = -A*X ; < ■ (96) 

Two of the wave functions (|9~5"D for i = 3 and 4 are therefore eigenfunctions of C^ N+2 ^ 
corresponding to the same eigenvalue —A. One obtains also the two other eigenvalues: 

n± = i[A(g + g- 1 ) ± y/A*(q-q-i)* + 4]. (97) 



Notice that if A = [m] q , then f2 = [m ± l] q , where we have used the notation fl3~3"|). The 
eigenfunctions corresponding to the eigenvalues Q ± are: 



^\u { ; i+1) + (n^-qA)4 l+1) ). (98) 



Using the eigenvalues and eigenfunctions ( |94] ) for N = 2 and the recurrence relations (|96] 



(P?|) and ( |9"8"D one can get all the eigenvalues and eigenfunctions of C^ N > for any even N. 



The eigenvalues are 



[n-m) q (m = 0, ... , JV) (99) 



with a degeneracy 



C#- m . (100) 



For q = 1 one recovers the spectrum given by ( see equation (|90|)). Notice that the 
similarity transformations used here are even matrices (see Sec. 3. 2), therefore the supertrace 
operation define by eq.(^) stays valid, 
b) N=2n+1 
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We start again with q — 1 and from the definition of C*( 2n+1 ) we have: 



2/i 



C (2n+1) = C (2n) + ± JJ ^. 



fe=l 



Using the equations (|86"D and (^)-(|90|), we get the diagonal form of (7( 2n+1 ): 



^(2n+l) 



101) 



(102) 



which has obviously the spectrum 

AT 



— — 2m (m = 0, 1, 2, . . . , n) 



with a degeneracy: 



(103) 



/~t2m 



(104) 



We now consider the case q ^ 1 . Instead of ( |101|) one has 



1 



C (2n+1) = C (2n) + p ^ ^ 

S + S' 1 



(105) 



The recurrence relation fl82|) stays valid as well as (|96])-(|97|) and (^). As opposed to N = 2n 
where we use the recurrence relations starting with C^ 2 \ for A" = 2n + 1 we start with C®. 
In the basis w|^(i = 1, 2, 3, 4), the matrix is 



; s 2 (s + s- 1 ) 0^ 



cm 



S + S~ l 



V 



[s + s 





- 1 ) s~ 2 









-1 
-1 



having obviously the following eigenvalues (eigenfunctions): 



(y/qu\ + u 





1" 




.2. 




"1" 


; 


.2. 



/ -, , V"-l 



The spectrum of C ,( ^ 2n+1 ' ) is therefore: 



(106) 



(107) 



24 



N 

2m 

2 



(m = 0,1,..., n) (108) 



with a degeneracy: 



(109) 



From the spectra of the matrix C, which plays the role of a transfer matrix (see equation 
(|§D), one can derive the mass spectra (the inverse of the correlation lengths) using equations 
(H) and ([1081) . For N = 2n we have 



2 J<? 



M m = ln rJV L2J ^ (m= l,...,n-2) (110) 
and for A^ = 2n + 1 we obtain: 



2 m k 



M m = ln l2lq (m = l,...,n). (Ill) 

[f - 2m\ q 

Therefore the system is always massive. It is interesting to note that in the large N limit, 
for q = 1, one obtains (see equations (|110|) -( |iTTf )): 



N N 
lim — M m = m (N = 2n); lim — M m = 2m. (N = 2n + 1). (112) 

N^oo 2 N^oo 2 

This implies that in the N — > oo, all correlation lengths diverge. Looking at the expressions 
( |112|) and having in mind that in conformal invariant theories one has similar expressions 
with iV substituted by L (N of O(N) replacing L, the size of the system, of the conformal 
invariant quantum chain) we would expect some similarity between both physics. The 
analogy, however is not so simple since the degeneracy of the level m also diverges (see 
equations (|100| ) and (|109|) ). An explicit calculation of the correlation functions in the large 
limit, which we didn't do, will clarify the issue. 

It is interesting to notice that for 0(3), spin S, (2S + l)-state quantum chain, VBS gives 
for the the smallest mass M 1; the following large S behavior Q: 

lims^M, = 1 (113) 

Comparing the equations ( |112[ ) with ( |113| ) we learn that in the asymptotic cases, the largest 
correlation length is given essentially by the number of states of the chain. 
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VI. CONCLUSIONS 



We have considered g-deformed O(N) symmetric, iV-state quantum chains defined by 
Hamiltonians given by equation. ([!]), ( |32| ) and (44). The symmetry is unbroken for free 
boundary conditions. For q ^ 1 the quantum group symmetry is broken for periodic bound- 
ary conditions. For q — 1, no symmetry might be left because of boundary terms which can 
be chosen as described in Sec. 4. Using algebraic methods, the ground-state wave functions 
for these quantum chains are known exactly for periodic, free and non-diagonal boundary 
conditions, they all correspond to energy zero. The wave functions are obtained using q- 
deformed Clifford algebras. These generalizes the construction of Affleck et al ||. Using 
the trace and supertrace operation in an auxiliary space, for N even and periodic boundary 
conditions, one obtains two ground-states one for momentum zero and one for momentum 
7i. This implies that even for a finite number of sites and periodic boundary conditions, the 
ground state is degenerate. For iV odd one obtains only translationally invariant ground- 
states. For free boundary conditions the degeneracy of the ground-state is 2 N ~~ 1 . This 
degeneracy is lifted by boundary terms. We have shown how to compute correlation func- 
tions and have derived all the correlation lengths. They are finite and diverge only for q = 1 
and L — > oo. What is the physical relevance of our results? For N = 4 we have shown in 



the Appendix A how the chain can be mapped into the extended Hubbard model [IT]. For 
all values of N one can map our quantum chains for obvious reasons into various ladder 
models [|l^] writing the on-rung interaction as a two-site interaction. If what one obtains is 
physically interesting remains to be seen. On the other hand the wave functions we obtain 
can be used as trial ground-state for more realistic models || . Can the procedure described 
here be extended to other quantum chains? The answer is yes. One can consider g-deformed 



Sp(N) symmetric chains. In this case instead of the Clifford algebra one gets [17| the q- 
deformed Heisenberg algebra as a tool to compute the wave functions. One can go even one 
step further and take quantum chains with the superalgebra Osp(M/N) as symmetry. In 
this case [|T7|] the algebra one uses to construct the wave functions is a combination of the 
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Clifford and Heisenberg algebras. These extensions are straightforward. Again, it is an open 
question if these extensions are interesting from a physical point of view. Last but not least, 
very simple quadratic algebras were discussed above, if more interesting ones (with X a and 
Y a in equations (|6|) and (0) unequal to zero) find their use in equilibrium problems, remains 
to be seen. They do in non-equilibrium problems. 
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APPENDIX A. FERMIONIC FORMULATION OF THE 0(4) QUANTUM CHAIN 

In this appendix we are going to present explicitly the Hamiltonian that corresponds to 



the N = 4 case. From (|32l) and ( 34]) the Hamiltonian is 



H = J2H k , (A.l) 

k 



4 

tt _ p(+) _ \ " -pa/3 77,7a jp5p _ 
U k - r h - 2^ L -yS- tj k %+l - 
a, /3, 7, (5=1 

^zt{(? + q-'M'EH, + EfEll, + EfEf +1 

,7-44 77.44 I 1 „[pHp22 _,_ 771 11 77133 , 7-t22 77'44 1 7Ti33 7Ti44 1 1 -If 77133 77.11 , 7714477133 
+^k %+lJ + Q\Pk h k+l + h k h k+l + h k h k+l + h k h k+l\ + Q Y^k h k+l + h k h k+l 

, 77144 77122 ,77-2277111 1 1 r 77121 77112 ,7711277121 1 7731 7713 1 7-ii3 77131 , 7742 77-24 

+E k E k+1 + E k E k+1 \ + [E k E k+1 + E k E k+1 + E k E k+1 + E k E k+l + E k E k+1 

, 77-24 77-42 , 77-4377-34 , 77-3477-43 l , n r 77122 77i33 , 7713377122 , 77 14 7741 1 77-2377132 
+E k E k+1 + E k E k+1 + E k E k+l \ + a 3 [E k E k+1 + E k E k+1 + E k E k+1 + E k E k+1 

+EfEf +l + E^ElU + a lE ?Et* +1 + a^E^ + a 2 [EfEi% x + EfE^ 

+El 2 E^ +1 + El^El 2 +l } + a 4 [EfEl 3 +1 + EfE]? +l + E^E^ + E^E^]}, (A. 2) 
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where 

q 3 q 2 q 1 q' 1 { . . 

«1 = ; | 2 ' " 2 = ~ l i 2 ' " 3 = 1 I 2 ' " 4 = ~ 1 i 2 ' "5 = : . 2 - (A.3) 

1 + 1 + <r 1 + 1 + 1 + q l 

It is also interesting to rewrite ( |A.2| ) in terms of spin-| creation and annihilation fermion 
operators on the lattice. This is done by making the following correspondence between the 



basis | a >j, a = 1, 2, 3, 4, in ( A.2Q , at each lattice point j, and the Fock representation 



|1 >j <"> 10 >i= I • • >i, |2 >;<-> c+ + |0 >,= | T ■ 

I 3 >J ^ 4-1° >i= I' i>i> I 4 4+4-1° >i= I t i>i • (A.4) 



Using this fermionic basis the Hamiltonian density ( A. 2 ) is given by 



H k = 1 _ 1 { ( c ta c k+i,a + h.c.)(l + t al n k - a + t a2 n k +i-a + t' a n k _ a n k+1 _ a ) 

q ~i~ q <r=+,- 

+J(5 fc ■ - n k n k+1 /4:) + t p (c£ >+ c£_c k+1>+ c k+1 _ + h.c.) + (g + g _1 ) - gn fc 
-g _1 ^fc+i + U\n k& n k „ + C/ r .n Jfe+1(+ njfc + i _ + ^ Kr,<r /7l v n fc+V 

ct,ct'=+,- 

+ [Va 'n kt+ n k+1 -n k+1<+ + K 3 V y n fe ,_n fc+ i _n fc+ i i+ + V3 'n k ^n kt+ n k+1 , + 

+V^ ) n k _n fci+ n fc+ i _] + ) _n fc)+ n fc+ i _n fc+ i,+}, (A. 5) 



where 



9+2 1 + 2g2 / _ / _ o 2g 

*-i - *+i - ~^^T' *-a - *+2 - - : + 2 ' *--*+- 3, J -2t p - Y^f' 

Ui = -^, U r = -^, V ++ = V^ = q + q-\ V + _ = V_ 2,1 



1 + g 2 1 + g 2 ^ 3 ' 3 ' T ~ '" r 1 + 

Vi 1) = Vi*> = -q-\ V 3 ^ = V^ = -q, Vt = q + q~\ (A.6) 



In ( |A.5| ) appear the density operators n fc (T = c~l a c k ^ and = + + n fc _ at the site k. The 



magnetic spin-spin interaction (coupling J) in ( A.5 ) is derived from the relation 



Y c k ,a c k+i,a' c k,a'C k+ i,a = 2(S k • Sfc+i - n fc ra fc+ i/4) + n k)+ n k+ i _ + rife _nfc+i )+ , (A. 7) 
where = |(7fe, and er = (cr", cr y , cr*) are the spin-| Pauli matrices. 
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The Hamiltonian ( |A.5|) belongs to the class of extended Hubbard models considered 



in the recent literature JTTJ. Beyond the magnetic interaction (coupling J) we also have 
non-diagonal interactions that correspond to single particle correlated hopping (couplings 
^o-i, tr2, a = as we U as P a i r hopping terms (coupling t p ). The static interactions 
are given by the diagonal terms. The couplings Ui and U r give us the on-site Coulomb 
interaction, and the interactions V a>a t (a, a' = ±), V^ a \ (a = 1, ... ,4) and V4 give us the 
two- three- and four-body static interactions, respectively. 

We should notice that the Hamiltonian ( |A.5|) conserves separately the total number of up 
spins n + and down spins n_. Consequently for free boundary consitions we may construct, 
using the the algebraic method, zero-energy eigenfunctions *f? n+)n _, for each sector labelled 
by n + and n_ (n + , n_ = 0, 1, . . . , L), i. e., 

L 



^ —V 

^n_|_,n_ ' n_)_,n_ 



J| (x! + x 2 4 t+ + x 3 4_ + x 4 cj + cj_) |0 >, 



(A.8) 



where P n+ ,n_ projects out states which do not have n + spins a = + and n_ spins a = — 
(see equation (|54"D). 

APPENDIX B. CORRELATION FUNCTIONS FOR PARITY VIOLATING 

OPERATORS (N EVEN) 

We would like to show how to compute the correlation functions 

< 0,-|P r Q s |0,+ > 

4r,s = ^ (B.lj 

where 

Z-,+ =< 0,-|0,+ > (B.2) 

which appear for and L even and periodic boundary conditions when the vacuum is 
degenerate. Here 

|0, + >= Tr(x Ql . . . x aL )u ai ...u aL (B.3) 
29 



corresponds to the parity +, momentum zero wave function, 

< 0, -| = Str(y 01 . . . y h )u T Pl . . . u t Pl = Tr( Jy Pl . . . yp L )u T Pl . . . ttJ L (B.4) 

corresponds to the parity — , momentum 7r wave function. The matrix J is defined in 
equation (fj5|). The action of the operators P and Q is shown equation (^). For obvious 
reasons, one of the two operators P or Q has to break parity. It is easy to show, using the 
definitions given by equations (|23|) and (P5|) that we have 



< 0, -|P r Q s |0, + >= Ti(DC r ~ 2 PC s - r ^QC L - s ) (B.5) 



and 



< 0,-|0,+ >= Tr(DC L - 1 ) (B.6) 



where 



v 



d = J2 x « ® J y*- ( R7 ) 

a=l 

Obviously the correlation lengths appearing for this type of correlation functions are the same 
as for the parity conserving operators where one computes quantities like <0,+|...|0,+ > 
or < 0, -I . . . 10, - >. 
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